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Resonance fluorescence arises from the inter¬ 
action of an optical fleld with a two-level sys¬ 
tem and has played a fundamental role in the 
development of quantum optics and its applica¬ 
tions. Despite its conceptual simplicity it entails 
a wide range of intriguing phenomena, such as 
the Mollow-triplet emission spectrum[l], photon 
antibunching[2] and coherent photon emission[3]. 
One fundamental aspect of resonance fluores¬ 
cence, reduced quantum fluctuations in the single 
photon stream from an atom in free space, was 
predicted more than 30 years ago[4]. However, 
the requirement to operate in the weak excitation 
regime, together with the combination of modest 
oscillator strength of atoms and low collection ef- 
flciencies, has continued to cast stringent experi¬ 
mental conditions for the observation of squeezing 
with atoms. Attempts to circumvent these issues 
had to sacriflce antibunching due to either stim¬ 
ulated forward scattering from atomic ensembles 
[5, 6] or multi-photon transitions inside optical 
cavities [7, 8]. Here, we use an artiflcial atom 
with a large optical dipole enabling 100-fold im¬ 
provement of the photon detection rate over the 
natural atom counterpart [9] and reach the neces¬ 
sary conditions for the observation of quadrature 
squeezing in single resonance-fluorescence pho¬ 
tons. Implementing phase-dependent homodyne 
intensity-correlation detection[9 12], we demon¬ 
strate that the electric fleld quadrature variance 
of resonance fluorescence is 3% below the fun¬ 
damental limit set by vacuum fluctuations, while 
the photon statistics remain antibunched. The 
presence of squeezing and antibunching simul¬ 
taneously is a fully nonclassical outcome of the 
wave-particle duality of photons. 

The minimum fluctuations in any quantum measurement 
of canonically conjugate variables such as position and 
momentum are bound by the Heisenberg uncertainty 
principle. Although this principle cannot be violated, 
the fluctuations of a single variable can be reduced below 
this minimum value at the expense of enhancing the 
fluctuations of the conjugate variable. The most widely 
explored realization of this nonclassical phenomenon is 
squeezed light [13] where the quadrature operators Xi 
and X 2 of the electric field are the canonically conjugate 
operators. Relying inherently on the quadratic depen¬ 
dence on the bosonic creation and annihilation operators 
in the Hamiltonian, squeezed light is typically generated 


using intense lasers and macroscopic nonlinear-optical 
media [14]. This form of squeezed light has multiple 
applications in the field of quantum optics[15], one 
prominent example being interferometry with reduced 
quantum noise[16, 17]. 

In 1981, Walls and Zoller predicted that the quadratic 
form of the Hamiltonian is not a requirement and that 
quadrature squeezing can also be obtained via a radi¬ 
cally different approach: the interaction of a two-level 
system with a resonant light field, as described by the 
Jaynes-Cummings Hamiltonian [4]. The fluctuations 
in one quadrature, quantified by their variance, can be 
reduced up to a theoretical maximum of 12.5% lower 
than vacuum fluctuations, while the intensity statistics 
remain antibunched. Unlike its nonlinear-optics coun¬ 
terpart this unique form of squeezed light stems from a 
buildup of atomic coherence which, once mapped onto 
the emitted field, results in the creation of coherences 
between the n = 0 and n = 1 Fock states in the weak 
excitation regime. Higher number states are excluded by 
photon antibunching or equivalently, by the fermionic 
nature of atomic operators. The simultaneous presence 
of antibunching and squeezing is an intriguing yet 
counter-intuitive effect, since single photons do not 
have a well-defined phase. It is the coherence with the 
zero-photon, i.e. vacuum, component that allows for 
phase-dependent effects such as the squeezing discussed 
here. 

The two-level system we use in this work is a voltage- 
controlled semiconductor quantum dot (QD)[18] 
positioned under a solid immersion lens for enhanced 
photon collection efhciency [19]. Typically allowing 
photon detection rates exceeding well above a million 
photons per second, these artificial atoms obviate the 
immediate need for cavity coupling and consequently 
allow for the experimental realisation of an isolated, 
weakly excited two-level system treated in ref. [4]. With 
large oscillator strength, high internal quantum efficiency 
and negligible decoherence, semiconductor QDs enabled 
recent observations of key phenomena in quantum optics, 
such as antibunching [20, 21], formation of dressed states 
[22-24], generation of entangled photon pairs [25-27] 
and, particularly relevant for this work, coherent single 
photon generation via weak excitation [28, 29]. The 
strong transition dipole moment of the QD compared to 
single atoms is also the key enabler in the detection of 
squeezing in resonance fluorescence, since signal to noise 
ratio is fundamentally tied to the photon detection rate. 

To generate resonance fluorescence, we excite the 
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FIG. 1. Homodyne intensity-correlation detection, a, Schematic illustration of the homodyne intensity-correlation setup. 
LP: linear polariser, PBS: polarising beam splitter, ND: neutral density filter, A/2: half-wave plate, b, Intensity of the SL field 
on a single detector as a function of the interferometer phase at an excitation power of s = 0.1. At this excitation power, each 
detector of the HBT setup records 1.6 • 10® events/s from resonance fluorescence contribution alone, c, Dipole phase offset in 
interference pattern as a function of detuning between laser and QD frequency, d, e, Intensity autocorrelation measurement 
with LO path blocked (d) and unblocked (e). In the blocked case, the measurement gives the expected antibunching, 

evidence of a single two-level system, regardless of (j>. The ordinate in panels d and e shows the coincidences in units of 
count rates for comparison. For unblocked LO (panel e), the interference between LO and QD fields leads to phase-dependent 
correlations, some of which contain the quadrature variance of the QD field. 


TT■'"-polarised neutral exciton transition of a single 
QD resonantly at 970 nm using a frequency-stabilised 
tunable laser (Fig. la). Resonance fluorescence and the 
reflected laser are separated by a polarising beam splitter 
and, after imparting a relative phase through a path 
length difference A£, recombined via a non-polarising 
beam splitter. One of the outputs of this beam splitter 
contains the superimposed light (SL) field 

where the subscript LO (RF) indicates local oscillator 
(resonance fluorescence) and the relative amplitude and 
phase (j) of the LO and RF fields can be tuned inde¬ 
pendently in the experimental scheme illustrated in Fig. 
la. Reflection and transmission coefficients and all other 
relative phases due to the optical setup are included in 
the field amplitudes and the phase (j). Figure lb shows 
the intensity measured on a single detector as a func¬ 
tion of the interferometer-induced phase kAi, where k 
is the wavenumber of the fields. The phase due to the 
dipolar response of the transition, which is determined 
by the relative detuning between the excitation laser and 


the transition frequency, is also included in </); Fig. Ic 
displays the measurement of the detuning dependence of 
this additional phase. 

The amplitude of a light field, where (j) is the 

relative phase with respect to a coherent reference field, 
can be represented in the phase space of the conjugate 
variables, Xx and X 2 , via E{(j)) ex X{(j)) = {Xicoscj) + 
X 2 sin (j)). These quadratures are the analogs of the di¬ 
mensionless position and momentum operators and their 
variances, quantifying the quantum fluctuations of the 
electric field, are subject to a similar uncertainty rela¬ 
tion: 

AXfAX^>l/16. (2) 

To demonstrate a squeezed quadrature, i.e. AX{(j))‘^ < 
1/4, we implement the experimental setup proposed in 
Ref. 12, which provides a direct and convenient link be¬ 
tween time-correlated two-photon detection and the vari¬ 
ances of field quadratures. To detect the variances of the 
RF field quadratures Ai_ 2 , we perform an intensity auto¬ 
correlation on the SL field ifsL using a Hanbury Brown 
and Twiss correlation setup (HBT). The unnormalised 


2 















0.15 



FIG. 2. Phase-dependent quadrature variance of resonance fluorescence, a, In-phase (orange) and out-of-phase 
(blue) normally ordered autocorrelations of the electric field quadrature fluctuations (AX(<(), 0)AX(()), r)): ^ for high and 

low power excitation conditions are displayed in top and bottom panels, respectively. Negative values in panel a (bottom) 
verify squeezing of the in-phase electric field variance, b, Full phase dependence of the quadrature variances (zero-time delay 
of the autocorrelations), for different excitation powers. A measurement of coherent laser quadratures provides a reference for 
the vacuum limit of zero (green circles, rightmost panel). Solid curves in all panels are theoretical simulations using a two-level 
Master equation for the corresponding experimental conditions. 


second-order correlation function of the SL field. 


AX{cl>,t): 


G[l,(i,i+r) = (4-)(t)^(-)(t + r)i;«(t + r)s(;)(t)). 

( 3 ) 

produces the well-known antibunched second-order corre¬ 
lation function of the resonance fluorescence field, 
in the absence of a local oscillator, showing the single 
photon nature of resonance fluorescence [2]. The solid 
red circles in Fig. Id display this behaviour for an ex¬ 
citation power of s = P/Psut = 0.1, where the satura¬ 
tion power Psat yields half of the maximum attainable 
resonance fluorescence intensity. The black curve is the 
theoretical prediction obtained with a two-level Master 
equation and includes the detector response function, as 
well as sub-linewidth spectral wandering of the QD tran¬ 
sition frequency [30, 31]. In the presence of the LO held, 
Gjotai displays a strong dependence on phase 0, as shown 
in Fig. le for (j) = 0, 7r/2 and tt. While the coincidence 
rate at long time delays changes with (j) by more than an 
order of magnitude, the correlation behaviour at short 
time delays evolves from a dip to a peak as a function of 
(f). 

( 2 ) 

The total correlation function contains hve terms 

with [pLor for n ranging from 0 to 4. We can separate 
their contributions via their unique dependence on the 
time delay (r) and relative phase (</>), as well as via di¬ 
rect measurement of the zeroth-order contribution (Fig. 
Id) (see SI). The second-order contribution is directly 
proportional to the normally ordered autocorrelation of 


AGf (r) cx ( : AX(<(>, 0) A1((/., t) : ) . (4) 


The zero-delay value 
ordered variance ( : 


of Eq. 4 hence yields the normally 
^AX{(j))^ : \ of the electric held 


quadrature [12]. This variance is zero for vacuum and for 
coherent helds and the existence of quadrature squeezing 
is manifested in a negative-valued normally ordered vari¬ 
ance: 


(aA(</.))': ^ <0. (5) 

Figure 2a presents the autocorrelation of the in-phase 
= 0) and in-quadrature {(j) = 7r/2) huctuations for 
high excitation power (s = 3, upper panel). The nor¬ 
mally ordered variance (r = 0) of resonance huorescence 
in the high power regime is positive valued regardless of 
the phase. This indicates that the quantum huctuations 
are enhanced above the vacuum level, as expected. In 
stark contrast, the low power regime (s = 0.1), shown in 
the lower panel, yields negative values for the normally 
ordered variance for (j) = 0. This reduction of quantum 
huctuations below the vacuum limit is the verihcation of 
quadrature squeezing in this measurement. As dictated 
by the Heisenberg uncertainty relations, this squeezing is 
accompanied by increased huctuations, i.e. antisqueez¬ 
ing, in the other quadrature. Both features decay on a 
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timescale of the order of the excited state lifetime [32]. 
Figure 2b shows how the normally ordered variance 
evolves as a function of (j), for excitation powers rang¬ 
ing from s = 0.05 (leftmost panel) to s = 3 (rightmost 
panel). A measurement with a weak laser of similar in¬ 
tensity is displayed in the rightmost plot as green circles. 
As expected for a coherent state, this measurement yields 
a normally ordered variance of zero, independent of </>. 
The squeezing vanishes [4] for s > 1, yielding larger fluc¬ 
tuations than vacuum for any cj). While we observe the 
phase dependence of the quadrature variance at all exci¬ 
tation powers, the window of opportunity for measuring 
negative values is restricted to a very small ^-range in 
the low-excitation regime, i.e. s < 1, highlighting the 
challenges for the observation of squeezing in resonance 
fluorescence since its prediction. 

Figure 3a summarises the power dependence of the nor¬ 
mally ordered quadrature variance extrema. The solid 
blue (red) curve represents the theoretically predicted 
behaviour of the in-phase (in-quadrature) field variance 
for an ideal two-level system. The maximum possible 
squeezing is limited to 12.5% (0.58 dB) below vacuum 
fluctuations at s = 0.36. The dashed curves depict how 
the ideal two-level system behaviour is modified due to 
the combined effects of finite timing resolution of our de¬ 
tection system and phase uncertainty of our interferome¬ 
ter (see Methods). All variance extrema we measure are 
commensurate with these predictions confirming that de¬ 
viation from the solid curve is predominantly of technical 
nature. The maximum degree of squeezing we measure is 
3.1 ± 1% (0.14 dB) below vacuum noise at an excitation 
power of s = 0.1. This value corresponds to 40% of the 
theoretically obtainable limit set by the blue solid curve 
at this excitation power. 

The transformation of the state of light with excitation 
power is best visualised by the calculated Wigner func¬ 
tions presented in Fig. 3b. The spread of these phase- 
space distributions along a given polar angle, </>, is indica¬ 
tive of the variance of the corresponding held quadrature 
X{(j)) = XiCOs{4>) + A 2 sin (</>). The Wigner function 
for the vacuum state (top left) shows a symmetric form 
with no phase dependence. At intermediate powers (top 
right and bottom left), the symmetry breaks down and 
a 0-dependence arises in the spread of the Wigner func¬ 
tion, linked to the generation of atomic coherence (see 
SI). This phase dependence, in combination with the an¬ 
tibunched nature of resonance huorescence, leads to a 
reduced variance of the electric held for a phase angle of 
0 = 0 (c.f. Wigner function at s = 0.36). In the high 
power regime (s —>■ 00 ), the held becomes a statistical 
mixture of n = 0 and n = 1 Fock states and the steady- 
state phase dependence is lost completely. 

We have shown that resonance huorescence from a two- 
level system can comprise a single photon stream with 
below-vacuum quantum huctuations of the held. While 
this appears counterintuitive owing to the impossibility 
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FIG. 3. Excitation-power dependence of quadrature 
squeezing, a, The measured (symbols) and the theoreti¬ 
cal (solid curves) normally ordered quadrature variances as 
a function of the excitation power s is shown for 0 = 0 and 
0 = 7r/2. The dashed curves illustrate the effect of phase noise 
and hnite timing resolution in our experiment on an otherwise 
ideal two-level system (see SI), b, Wigner functions for dif¬ 
ferent excitation powers. The transition from vacuum state 
(s = 0) to a mixture of vacuum and single photon Fock state 
(s —>■ 00 ) displays non-symmetric forms at intermediate exci¬ 
tation regimes (s = 0.36 and s = 10). The panel for s = 0.36 
displays clearly that the spread of the Wigner function along 
Xi is less than that of vacuum, a manifestation of quadrature 
squeezing. The dashed lines depict the contour at 50% of the 
maximum value at each power. 


of associating a well-defined phase to single photons, the 
probabilistic nature of coherent photon scattering in the 
weak excitation regime allows the emitted photons to be 
in a coherent superposition of Fock states |0) and |1). 
The emergence of phase correlations in this regime en¬ 
dows resonance fluorescence with the coexistence of pho¬ 
ton antibunching and quadrature squeezing. Our simul¬ 
taneous observation of these two phenomena can there- 
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fore be interpreted as a quantum mechanical manifesta¬ 
tion of the complementary particle and wave natures of 
light, respectively, with no classical analogs. 

We gratefully acknowledge financial support by the 
University of Cambridge, the European Research Council 
ERC Consolidator Grant Agreement No. 617985 and the 
EU-FP7 Marie Curie Initial Training Network S3NANO. 
C.M. gratefully acknowledges Clare College Cambridge 
for financial support through a Junior Research Fellow¬ 
ship. We thank E. Clarke, M. Hugues and the EPSRC 
National Centre for III-V Technologies for the wafer and 
C. Baune, R. Moghadas Nia, A. Ourjoumtsev, W. Vogel 
and H. J. Carmichael for fruitful discussions. 

Correspondence and requests for materials should be 
addressed to M.A. (ma424@cam.ac.uk). 

METHODS 

Interferometer: A frequency- and power-stabilised 
single-mode laser is used to resonantly excite the neutral 
exciton transition of the QD. The emitted photons are 
collected in a confocal dark-field microscope, where the 
laser is separated from the QD emission by means of two 
crossed polarisers. The second polariser is implemented 
as a polarising beam splitter (PBS), which enables the 
use of the attenuated laser field as local oscillator (LO). 
The light field in the QD arm of the interferometer con¬ 
sists typically of < 1% laser photons and > 99% RE pho¬ 
tons, we therefore neglect the laser background in the 
QD-photon mode. Likewise, any QD photons in the LO- 
output can be neglected because the excitation laser in¬ 
tensity before attenuation is several orders of magnitude 
larger than the RE intensity. In the fringe measurement 
in Fig. la the spatial path difference A£ « 11 cm is kept 
constant while laser and QD frequency are tuned con¬ 
tinuously to change the interferometer phase. This form 
of phase control is enabled by the tuneability of the QD 
resonance via the DC Stark effect [33] and increases the 
long-term stability of the interferometer, which contains 
no moving parts. The visibility of the interferometer for 
high power laser light amounts to near unity but is re¬ 
duced in Fig. lb to 73.8% due to incoherent photon 
emission as well as an inadvertent mismatch of RE and 
LO intensities. Additionally, the visibility is reduced for 
low count rates [9], which makes the use of a bright sin¬ 
gle photon source and high photon collection efficiencies 
crucial for our experiments. The collection efficiency of 
our optical setup is 1%. This value is calculated from 
the obtained count rates and the 0.58—ns excited state 
lifetime of the QD used in this work. 

Postselection: The intensities of RE and LO fields are 
kept equal in all of our measurements. To ensure the 
absence of laser photons in the QD mode, the laser back¬ 
ground is measured once a minute for 2 s. To this end, 
the LO path is blocked and the QD, which is embed¬ 


ded in a Schottky diode structure, is tuned oft resonance 
via the quantum confined Stark effect (QCSE) [33]. Fur¬ 
thermore, the intensity of the QD emission is monitored 
continuously during the measurements to detect spectral 
wandering of the QD transition [30, 34]. This is done by 
filtering out the phonon sideband (PSB) and detecting it 
on a third single photon detector [35]. The measured cor- 
relation histograms are saved once a minute. 

We perform a postselection of histograms with a thresh¬ 
old on the mean PSB count rate and another threshold 
on the measured laser leakage in the QD arm. In the ex¬ 
periments shown in Fig. 2, the laser is kept on resonance 
but the relative phase is not actively controlled in the 
interferometer. Instead, phase-dependent measurements 
are performed by using the individual detector intensi¬ 
ties as a measure of the interferometer phase and relying 
on the wandering of the phase due to temperature drifts 
on timescales of typically > 30 min/7r. In order to bin 
the data we perform a reference measurement of the in¬ 
terference fringes by scanning the laser frequency while 
keeping the QD on resonance using the QCSE. An exam¬ 
ple measurement of the interference fringes obtained in 
this way can be seen in Fig. lb. In order to have equal 
sized phase bins, we use intensity bins of varying size 
proportional to the derivative of a cos^((()/2) function. 
We note that this phase binning is sign-invariant, i.e. it 
cannot distinguish between positive or negative phases 
and therefore bins data into values between 0 and J-tt. 
This does not impact our measurement as all correlation 
functions are symmetric in phase around 0. The data 
points shown at negative phases in Fig. 2b are measured 
between (j) = 0.57r and (j) = tt, and have been shifted by 

— TT. 

Reduction of measured degree of squeezing: Al¬ 
though the conditioning nature of the measurement 
should render it robust against low detection efficien¬ 
cies, this is not true in practice, and several effects re¬ 
duce the measured degree of squeezing compared to the 
theoretical limit for resonance fluorescence. Low photon 
numbers in the interferometer reduce the fringe visibility, 
but this only affects the signal-to-noise ratio in our mea¬ 
surements. However, low count rates also lead to higher 
shot noise which can increase the error in the phase bin¬ 
ning protocol. This in turn leads to a decrease in the 
detected degree of squeezing by introducing mixing be¬ 
tween quadratures [36]. Other sources of phase noise in¬ 
clude spectral wandering of the QD transition leading to 
fluctuations in the dipole phase, and any interferomet¬ 
ric instability on timescales shorter than the histogram 
saving time. Finite timing resolution of the correlation 
setup also leads to a decreased visibility of the features at 
0 time delay and further reduce our measured value for 
squeezing. We have independently measured the timing 
resolution of the HBT setup with a mode-locked pulsed 
laser source (< 3 ps pulse width) for different mean count 
rates. The extent of phase-noise from different sources is 
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harder to quantify and is used as a fitting parameter in 
the theoretical curves in Fig. 2. 
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